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We present the first large-acceptance measurement of event- wise mean transverse momentum (pt) 
fluctuations for Au-Au collisions at nucleon-nucleon center-of-momentum collision energy ^/sjvjv 
= 130 GeV. The observed non-statistical (p t ) fluctuations substantially exceed in magnitude fluc- 
tuations expected from the finite number of particles produced in a typical collision. The r.m.s. 
fractional width excess of the event-wise (pt) distribution is 13.7 ± O.l(stat) ±1.3(syst)% relative to 
a statistical reference, for the 15% most-central collisions and for charged hadrons within pseudo- 
rapidity range \n\ < 1, 2n azimuth and 0.15 < pt < 2 GeV/c. The width excess varies smoothly 
but non-monotonically with collision centrality, and does not display rapid changes with centrality 
which might indicate the presence of critical fluctuations. The reported (p t ) fluctuation excess is 
qualitatively larger than those observed at lower energies and differs markedly from theoretical ex- 
pectations. Contributions to (p t ) fluctuations from semi-hard parton scattering in the initial state 
and dissipation in the bulk colored medium are discussed. 



PACS numbers: 24.60.Ky, 25.75.Gz 



3 



I. INTRODUCTION 

Fluctuation analysis of relativistic heavy ion collisions 
has been advocated to search for critical phenomena near 
the predicted hadron-partonphase boundary of quantum 
chromodynamics (QCD) 0,0,0. Nonstatistical fluctua- 
tions (excess variance beyond statistical fluctuations due 
to finite particle number), varying rapidly with collision 
energy, projectile size or collision centrality and inter- 
preted as critical fluctuations, could indicate a transition 
to a quark-gluon plasma 0, 0, L| • Nonstatistical fluctu- 
ations could also appear in systems incompletely equili- 
brated following initial-state multiple scattering (Cronin 
effect Q and minimum-bias hard parton scattering - 
minijets or as an aspect of fragmentation of color 
strings produced in nucleon-nucleon collisions or the bulk 
medium in A- A collisions. Study of nonstatistical fluctu- 
ations and the correlations which produce them are cen- 
tral aspects of the Relativistic Heavy Ion Collider (RHIC) 
research program. The specific goal of the present work 
is to determine the magnitude and collision centrality 
dependence of nonstatistical fluctuations in momentum 
space at large momentum scales using the largest angular 
acceptance detector available at RHIC. 

The dynamical representation of relativistic nuclear 
collisions can be separated into transverse (perpendicular 
to the beam axis) and longitudinal (parallel to the beam 
axis) phase spaces. In this paper we focus on transverse 
phase space, including transverse momentum magnitude 
p t and momentum azimuth angle <p, within relatively 
small pseudorapidity rj intervals. Assuming rapid lon- 
gitudinal (Bjorken) expansion of the collision system Q, 
separate r\ intervals can be treated as quasi-independent 
(causally disconnected) dynamical systems. In this anal- 
ysis we calculate the event-wise mean transverse momen- 
tum for each collision event within a detector kinematic 
acceptance 

I N 

fa) = Jf^Pt.i, (!) 

i=l 

where i is a particle index and N represents the measured 
charged-particle multiplicity within the detector accep- 
tance for a given collision event. Quantity (p t ) is mono- 
tonically related to the 'temperature' of the event-wise pt 
distribution, plus any collective transverse velocity of the 
collision system. The distribution of (pt) over a collision 
event ensemble, especially any excess variance of this dis- 
tribution beyond what is expected for purely statistical 
fluctuations, reflects the underlying dynamics and degree 
of equilibration of heavy ion collisions. 

Some aspects of heavy ion collisions produce correla- 
tions/fluctuations which depend on the relative charge of 
a charged hadron pair B, 0, 01 , including quantum and 
Coulomb correlations [l(J, resonance decays, color-string 
fragmentation (e.g., charge ordering along the string 
axis 0, 0|), and minijet fragmentation. Charge de- 
pendent combinations for pion pairs can be directly re- 



lated to isospin components. For non-identified charged 
hadron pairs in the collisions studied here, which are 
dominated by pions but include other charged hadrons 
{e.g. protons, kaons and their antiparticles), the re- 
lation to isospin remains useful but becomes approxi- 
mate. To isolate the different isospin aspects of fluctua- 
tions and correlations we measure separately the like-sign 
(LS) and unlike-sign (US) charge-pair contributions and 
also form charge-independent (CI) and charge-dependent 
(CD) combinations, with CI = LS + US (approximately 
isoscalar) and CD = LS — US (approximately isovector) 
respectively. 

In this paper we report the first large-acceptance mea- 
surement of (pt) fluctuations at RHIC using the STAR 
detector. Results are presented for unidentified charged 
hadrons using I83k central and 205k minimum-trigger- 
bias ensembles of Au-Au collision events at ^/snn = 130 
GeV (CM energy per nucleon-nucleon pair). Experimen- 
tal details and the observed (pt) distribution for central 
events are presented in Sees. UTI and lllll Quantities used 
to measure nonstatistical (pt) fluctuations are discussed 
in Sec. II VI and the Appendix. Results and discussion are 
presented in Sees. W\- 1 Villi the observed large excess of 
(pt) fluctuations at RHIC is compared to other measure- 
ments and to theoretical models, including hard parton 
scattering in the initial state and/or hadronic rescatter- 
ing. Conclusions are presented in Sec. IIXI 



II. SUMMARY OF EXPERIMENT 

Data for this analysis were obtained with the STAR 
detector 0] employing a 0.25 T uniform magnetic field 
parallel to the beam axis. Event triggering with the 
central trigger barrel (CTB) scintillators and zero-degree 
calorimeters (ZDC) and charged-particle kinematic mea- 
surements with the time projection chamber (TPC) are 
described in H3|. TPC tracking efficiency was deter- 
mined to be 80 - 95% within \f]\<l and p t > 200 MeV/c 
by embedding simulated tracks in real-data events |14| . 
and was uniform in azimuth to 3% (r.m.s.) over 2ti. Split- 
track removal required the fraction of valid space points 
used in a track fit relative to the maximum number pos- 
sible to be > 50%. A primary event vertex within 75 cm 
of the axial center of the TPC was required. Valid TPC 
tracks fell within the full detector acceptance, defined 
here by 0.15 < p t < 2.0 GeV/c, \t)\ < 1 and 2-7T in az- 
imuth. Primary tracks were defined as having a distance 
of closest approach less than 3 cm from the reconstructed 
primary vertex which included a large fraction of true 
primary hadrons plus approximately 7% background con- 
tamination |14|. 

Two data sets were analyzed: 1) 183k central triggered 
Au-Au collision events constituting the 15% most-central 
collisions as determined by scintillator hits in the STAR 
CTB and 2) 205k minimum-bias collision events triggered 
by ZDC coincidence. The latter events were divided into 
eight centrality classes based on TPC track multiplicity 
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in 1 77 1 < 0.5 0| , the eight event classes comprising ap- 
proximately equal fractions of the upper 87±2% of the 
Au-Au total hadronic cross section. 



III. MEAN p t DISTRIBUTION 

The frequency distribution of event-wise (pt) f° r 183k 
15% most-central collision events is first studied graph- 
ically. The data histogram is compared to a statisti- 
cal reference distribution and is examined for evidence 
of anomalous event classes which could indicate either 
novel collision dynamics pj or experimental anoma- 
lies. The event-wise (p t ) data distribution is shown 
as the histogram in the upper panel of Fig. ^ Those 
data, representing 80 ± 5% of the true primary parti- 
cles within the acceptance, were binned using quantity 
((pt) - Pt)/(&pt/Vn), where p t and cr| t are respectively 
the mean and variance of the inclusive pt distribution 
of all accepted particles in the event ensemble and n is 
the event-wise multiplicity within the defined acceptance. 
That choice of event-wise random variable rather than 
(pt) is explained as follows. 

For independent particle pt samples from a fixed par- 
ent distribution (no nonstatistical fluctuations) the r.m.s. 
width of the frequency distribution on (p t ) is itself de- 
pendent on event multiplicity n as o~p t / ' yfn (central limit 
theorem or CLT [lj, [lfj). The underlying purpose of 
this measurement is to determine an aspect of pt fluc- 
tuations which is independent of event multiplicity per 
se. If n is a random variable, a systematic dependence 
is introduced into the measured (pt) fluctuation excess 
through this CLT behavior of the width. To insure mul- 
tiplicity independence the basic statistical quantity must 
be formulated carefully. By normalizing the distribution 
variable with factor \/n/o~p t , the distribution width of 
the new variable is unity, independent of n, when fluc- 
tuations are purely statistical. The trivial broadening 
of the (pt) distribution for event ensembles with a finite 
range of event multiplicities is eliminated. The latter ef- 
fect can have significant consequences for relevant event 
ensembles (p-p, peripheral A-A and small detector ac- 
ceptance). This argument explains the variable choice 
for Fig. as well as the associated numerical analysis 
described in Sec. lIVI For the sake of brevity this normal- 
ized variable will in some cases still be referred to in the 
text as "(pt)". 

The precision of these data warrants construction of a 
statistical reference which accurately represents the ex- 
pected (pt) distribution in the absence of nonstatistical 
fluctuations. Because of its close connection to the cen- 
tral limit theorem (behavior under n-folding noted be- 
low) we can compactly and accurately represent the 
reference distribution with a gamma distribution 
We observe that the measured inclusive pt distribution 
is, for present purposes, well approximated by a gamma 
distribution with folding index ao = Ptl a § t ~ 2. Differ- 
ences between the gamma and inclusive pt distributions 




FIG. 1: Upper panel: Event frequency distribution on 
\/n((pt) — pt)/<Tp t (see text) for 80% of primary charged 
hadrons in ] -7-7 1 < 1 for 183k central events (histogram) com- 
pared to gamma reference (dashed curve), Monte Carlo ref- 
erence (solid curve underlying gamma reference), and broad- 
ened distribution (solid curve underlying data, not a fit - see 
text). Lower panel: Difference in upper panel between data 
and gamma reference (histogram) or between broadened dis- 
tribution and gamma reference (solid curve) normalized by 
the Poisson error \^N evt in each bin. 



in the higher cumulants due to p t acceptance cuts and 
physics correlations are strongly suppressed in the com- 
parison with the distribution in Fig.^by inverse powers 
of event multiplicity and are not significant for central 
Au-Au collisions. 

Because the n-folding of a gamma distribution is also 
a gamma distribution (representing an ensemble of inde- 
pendent n-samples of the parent gamma distribution or 
inclusive distribution) , the (p t ) reference distribution 
can be represented by [17| 



9n((Pt)) 



ao 
Pt 



-a n(p t )/pt 

T(a n) 



ao n- 



pt 



ao n—1 



(2) 



The corresponding gamma-distribution reference is indi- 
cated by the dashed curve in the upper panel of Fig. Q 
Parameter values used for this reference curve were deter- 
mined from the measured inclusive p t distribution as n = 
735 ± 0.2, p t = 535.32±0.05 MeV/c, and a Vt = 359.54± 
0.03 MeV/c, obtained from all accepted particles and not 
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corrected for p t acceptance cuts and inefficiencies. 

A reference can also be generated by a Monte Carlo 
procedure. An ensemble of n-sample reference events 
is generated with multiplicity distribution similar to the 
data. A reference event with multiplicity n drawn from 
that distribution is assembled by performing n random 
samples from a fixed parent pt distribution estimated by 
the interpolated inclusive pt histogram of all accepted 
particles from all events in the centrality bin. The result- 
ing Monte Carlo reference distribution is shown in Fig. 
upper panel by the solid curve underlying the dashed 
gamma reference curve. The agreement is excellent. The 
broadened distribution (solid curve) underlying the data 
in the upper panel of Fig. ^ is discussed in Sec. E] All 
curves are normalized to match the data near the peak 
value, emphasizing the width comparison which is the 
main issue of this paper. We observe a substantial width 
excess in the data relative to the statistical reference. 

Fig. n lower panel shows the difference between data 
and gamma reference normalized to Poisson standard de- 
viations in each bin, emphasizing the large statistical sig- 
nificance of the width excess. We observe no significant 
deviations (bumps) from the broadened distribution in 
Fig. ^ which might indicate anomalous event classes as 
expected in some phase-transition scenarios pj . It is also 
important to note that the entire event ensemble con- 
tributes to the width increase relative to the statistical 
reference, i.e., the excess width is not dominated by a 
subset of problematic events. We note that the distribu- 
tion in Fig.^cannot be corrected for background contam- 
ination and tracking inefficiency. The numerical analysis 
described in the next section allows such corrections. 



IV. MEASURES OF NONSTATISTICAL (p t ) 
FLUCTUATIONS 

Consistent with the argument presented above about 
eliminating dependence of fluctuation measures on multi- 
plicity variations within a centrality bin, we characterize 
the magnitude of nonstatistical (pt) fluctuations by com- 
paring the variance of distribution quantity y/n((p t ) —pt) 
from Fig. ^ to the variance cr| t of its reference distribu- 
tion. The difference between these two variances is rep- 
resented by 

Ao%™> = - [{pth -Pt} 2 - 4 t (3) 

i=i 

= 2a Pt A^-l (4) 

where e is the number of events in a centrality bin, j is 
the event index, rij is the number of accepted particles 
in event j, (ptjj is the mean pt of accepted particles in 
event j. Subscript pt : n emphasizes that this quantity 
measures variance excess due to fluctuations of pt rela- 
tive to event-wise fluctuations in multiplicity n (i.e., is 
not significantly affected by fluctuations in n itself) . Su- 
perscript (CI) indicates a charge-independent sum over 



all particles. Difference factor Aa Pt -.n defined in Eq. (0J 
is approximately equal to (p f ) fluctuation measure <& pt 
introduced previously 0, ^3 • 

Two issues motivate the definition of fluctuation mea- 
sure A<7pf^? in Eq. 10: I) (p t ) is the ratio of two ran- 
dom variables - a scalar pt sum and a multiplicity. Fluc- 
tuations in either variable contribute to fluctuations in 
the ratio. For an uncorrelated system with fluctuating 
multiplicity, ratio fluctuations go as 1/y/n, producing 
an apparent nonstatistical contribution to ratio fluctu- 
ation measures which are aimed at determining p t fluc- 
tuations. 2) Measures of nonstatistical fluctuations typ- 
ically involve (at least implicitly) a difference between 
variances evaluated at two different scales, where 'scale' 
in the present context refers to histogram bin sizes (e.g., 
on r] and <f>). Bins on 77 and </> are denoted respectively 
by 5rj and 5<p or generically by 5x. The detector accep- 
tance can define one scale, as in this analysis. The other 
relevant scale, both for the simulated events presented 
in the preceding section and in the variance measure- 
ments presented in Sec. [V] is the single particle scale in 
which the bins are always made small enough such that 
occupied bins contain a maximum of one particle. In 
general the scale is independent of the acceptance where 
scale < acceptance. The case of variance calculations for 
arbitrary scale is treated in the Appendix. Scale depen- 
dence of variance excess provides important information 
on the underlying two-particle correlations and is an es- 
sential feature of any nonstatistical fluctuation measure- 
ment such as those presented here, although the impor- 
tance of this point has not been fully appreciated in this 
heavy ion context. 

In the Appendix we show that the scale invariance 
of total variance, an expression of the CLT, motivates 
the quantity in Eq. ©. A<Tp£n (Sx) measures changes 
in variance stemming from two-particle correlations with 
characteristic lengths less than binning scale Sx ^(|. As 
a function of binning scale A<jpf :n (fa) is not dependent 
on an acceptance size (knowledge of its scale dependence 
may of course be limited by a finite detector acceptance) 
but can depend on absolute position of the acceptance in 
momentum space. 

Given the definition of <& Pt and Eq. J3J), Aop^n — 
(% t + °p t ? - Oft, and % t ~ Aa£% M Differ- 
ence factor Acrp^n and <& Pt are therefore comparable be- 
tween different analyses. Fluctuation measure a^ t dyn = 

((Pt,i - Pt)(Pt,j ~ Pt))^ H3 (bar denotes event average) 
is related to Ao^S? by a 2 pt dyn ~ Ao%$ /(N - 1) (N is 
the mean multiplicity) for approximately constant event- 
wise multiplicities. $ Pt and er^ d may include signifi- 
cant dependence on multiplicity fluctuations in the case 
of small bin multiplicities (e.g., for any bins within p-p or 
peripheral A-A events or for small-scale bins within cen- 
tral A-A events). Variance difference Aap^ minimizes 
this dependence compared to the preceding quantities. 
In Eqs. Q and the Appendix the summations over 
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particles have ignored charge sign. Aapl^ is a charge- 
independent (approximately isoscalar) quantity. By sep- 
arating contributions to Eq. J3J) into sums over (+) and 

(— ) charges, charge-dependent (CD) quantity AopfS, 
can be defined which measures the difference between 
contributions to (pt) fluctuations from like-sign pairs and 
unlikc-sign pairs. Using explicit charge-sign notation 
quantities Aa^n 



and Act 



2(CD) 



Pf-n 



are defined by 

iV(Ax)Aa^) = N(Ax) + Aa 2 pt ., n>++ 
+ 7V(A.x)_A^ t:ni _. 

+ 2JN{Ax) + N(Ax)-Aa. 



(5) 



2 

Pt:n,+- 



N{Ax)Aa^ 



N(Ax) + Aa 2 



p t :n,++ 
2 



(6) 



•V(A,-j Aaj„ : , 

2JN(Ax) + N(Ax).Aa 2 pt .^ + _ 



where N(Ax)± are the mean multiplicities for ± charges 
in acceptance Ax and N(Ax) is the mean total multi- 
plicity in Ax. Individual terms in Eqs. (JSJ and © are 
defined by 



A _2 

^"pt-n.ab 



n a {{pt)a - Pta) y/nb({pt)b ~ Ptb) 
1 : (7) 



' a 3t,a.^ab, 



where subscripts a and b represent charge sign, ab — 

++, ,+— or — h, the bar denotes an average over 

events, and 5 a b is a Kronecker delta. Difference factors 
A(Tp^ : n and Acr^^t) (approximately isoscalar and isovec- 
tor respectively) reported in the following sections are 
defined by 



Acr 2(CI) - 2cr- Acr (CI) 
A 2(CD) _ o a (CD) 



(8) 
(9) 



reference gamma distribution in Eq. (0) to the power 

er| t /(er| t + Aapl^P). This procedure, which would be ex- 
act for a gaussian distribution, increases the variance of 
the modified gamma distribution to the numerical value 
obtained from the data, preserves the mean, and agrees 
well with the relative peak heights of the data in the lower 
half of Fig. ^ The comparison in Fig. ^ then demon- 
strates that Aa Pt -n provides an excellent description of 
the event-wise (pt) distribution and its fluctuation excess. 
The corresponding r.m.s. width increase relative to the 
reference is 13.7 ± 0.1(stat)±1.3(syst)%. When extrap- 
olated to 100% of primary hadrons and no backgrounds 
A(jp^ : I r J CD ' > was estimated to be a factor 1.26 larger in 
magnitude for the 15% most-central events, resulting in 
a corrected charge- independent r.m.s. width increase of 
17±2(syst)%. 
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FIG. 2: Mean-pt difference factors Aa 



(CI) 



RESULTS 



and A<Tp^:n for 
205k minimum-bias Au-Au events at y^ijvjv = 130 GeV ver- 
sus relative multiplicity N/No Il4jl . which is approximately 
N V art/Nvart.max, the relative fraction of participant nucle- 
ons |2]J. Charge-independent (solid triangular points) and 
charge-dependent (open triangular points - multiplied by 3 
for clarity) difference factors include statistical errors only 
(smaller than symbols). Parametrizations (dashed curves), 
extrapolation of parametrizations to true primary particle 
number (solid curves), and systematic uncertainties (bands) 
are discussed in the text. Difference factors for the 15% most- 
central collision events are shown by the solid circle and open 
circle symbols. 



We apply Eqs. iJSJ - 10 to central collisions and to a 
minimum-bias ensemble. In all cases charge symmetry 

Aa 2 t . n++ ~ A(Tp t . n is observed within errors. For 

the 15% most-central events and full acceptance we ob- 
tain difference factors Aa p ^}n = 52.6 ± 0.3 (stat) MeV/c 
and Aap^i — —6.6 ±0.6 (stat) MeV/c (respectively the 
solid and open circular data symbols in Fig. [2J . Charge- 
independent values of <E> Pt and u ptdyn for the same data 
are respectively 52.6 ± 0.3 (stat) MeV/c and 52.3 ± 0.3 
(stat) (MeV/c) 2 (note units). Dependence on multiplic- 
ity fluctuations is negligible for this full-acceptance 15% 
most-central collision ensemble. 

The experimental value Aa p ^n = 52.6 MeV/c was 
used to determine the solid curves underlying the data 
histogram in the two panels of Fig. ^ by raising the 



Difference factors were also determined for eight cen- 
trality classes defined for the 205k minimum-bias events 
described in Sec. |H] Measured values of Aa p ^}n and 
Aap^n are shown in Fig. [2] by the upper and lower 
set of data symbols for CI and CD respectively, plot- 
ted for each centrality class vs its mean multiplicity N in 
\v\ < 0.5 (Sec.[nJ| relative to iVn, the minimum-bias mul- 
tiplicity distribution endpoint [2l| where N — 520 ± 5. 
Data are listed in Table U Plotted points, including 
statistical errors only (typically ±0.5 MeV/c), were fit- 
ted with parametrizations (dashed curves) which were 
then extrapolated by amounts varying from 1.17 to 1.26 
(for peripheral to central events respectively) to produce 
estimates for 100% of primary charged hadrons (solid 
curves). Au^'i has a very significant non-monotonic 
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dependence on centrality, but with no sharp structure. 
A<7p^n is significantly negative and approximately inde- 
pendent of centrality. <J> Pt and cr^^ (N — l)/2ap t agree 

with Acrp t: „ within statistical errors for the upper six 

(CD 

centrality classes, but both differ from Aerp t: „ and each 
other by much more than statistical uncertainty for the 
two most peripheral bins, as expected from their depen- 
dencies on multiplicity fluctuations. 

TABLE I: Centrality dependences of the measured charge in- 
dependent (CI) and charge dependent (CD) difference factors 
Acrp^n and Aop^„ plus the corresponding values extrapo- 
lated to 100% tracking efficiency. 



Centrality 


A<7 (01) 


(MeV/c) 


Aa iCb> 


(MeV/c) 


N/No 


<Tf„t \ / 


Np ar t b 


Data c 


Ext." 


Data c 


Ext." 


0.012 


87-76 


8.9 


22.8 


26.8 


-11.1 


-13.0 


0.033 


76-65 


19 


34.3 


40.4 


-6.9 


-8.1 


0.073 


65-54 


36 


39.6 


46.8 


-7.9 


-9.3 


0.14 


54-43 


64 


48.7 


57.6 


-7.4 


-8.8 


0.24 


43-33 


102 


51.3 


61.4 


-7.7 


-9.2 


0.38 


33-22 


153 


56.1 


68.0 


-7.0 


-8.5 


0.57 


22-11 


224 


54.4 


66.9 


-6.0 


-7.4 


0.84 


11-0 


320 


51.8 


65.1 


-6.6 


-8.3 



"Fraction of total hadronic inelastic cross section ranges in per- 
cent; values are ±2% uncertain 111 . 

'Estimates in |I4| were interpolated to centrality bins used here. 

statistical errors are typically ±0.5 MeV/c; systematic errors are 
±9%. 

^Difference factors extrapolated to 100% tracking efficiency and 
no secondary particle contamination. Uncertainties are ±12% and 
are dominated by systematics as discussed in the text. 

Systematic errors from uncertainty in two-track inef- 
ficiency, primary-vertex transverse position uncertainty, 
TPC drift speed/time-offset uncertainty, and conver- 
sion electron contamination were estimated by Monte 
Carlo [2^] as less than 4% of reported values. Stabil- 
ity of reported results against primary-vertex longitudi- 
nal position variation, momentum resolution, and TPC 
central membrane track crossing was determined to be 
5% of stated values. Systematic effects due to possi- 
ble time dependence in detector performance and effi- 
ciency were studied by analyzing sequential run blocks 
which were determined to be consistent within statistical 
error. Systematic error contributions due to azimuthal 
anisotropy in the event-wise primary particle distribu- 
tion (cos [2(0 — ^ r)] assumed where VP r is the event- 
wise reaction plane angle) combined with non-uniform 
azimuthal tracking efficiency was determined to be less 
than 1% of reported values using 0-dependent track cuts 
and measured efficiency maps. Nonprimary background 
(-7%) added ±7% systematic error due to uncer- 
tainty in its correlation content. Total systematic un- 
certainty for the A<Tp^ : „ and Aop^ : „ data in Fig. [21 and 
Table [I] is ±9%. Additional systematic error in extrapo- 
lation of AlTp^n and Aaffi to 100% of primary particles 
(±8%) is dominated by uncertainty in the actual primary 
particle yield |14| . Total uncertainty in extrapolated val- 



ues is about ±12% (shaded bands in Fig. 0). Systematic 
error in the most peripheral bin is larger by an addi- 
tional ~ ±1 MeV/c due to possible primary- vertex re- 
construction bias. Analyses of 30k central hijing Au-Au 
collision events both with and without STAR acceptance 
and event reconstruction effects yield consistent results 
for A<7 Pt ;n to within the statistical error (~10%) for these 
simulated events, which is well within our estimated sys- 
tematic error. 

Data in Fig. Hand Table [Q were not corrected for two- 
track inefficiencies, which would increase all results in a 
positive sense by up to 3 MeV/c. Variations (w 10%) in 
p t and tr| t with collision centrality were accommodated 
by independent analyses in small centrality bins. Monte 
Carlo |23 estimates indicate that combined corrections 
for quantum (HBT) and Coulomb correlations 01 j reso- 
nance (p°, u) decays, and p t centrality dependence (i.e., 
well known physical effects) would increase the absolute 
magnitudes of all data in Fig.|2land Table[I]by as much as 
« 6 MeV/c. Quantum and Coulomb correlations and res- 
onance decays originate in the final stage of the collision 
evolution and are not the main object of this study. Cor- 
recting AcTp^n for two-track inefficiencies plus the pre- 
ceding effects (not done for the data shown in Fig. and 
TableQ would cause the overall magnitude to increase by 

about 7 MeV/c. Similarly, corrections to Aa^' would 
cause it to become more negative by about 4 MeV/c. We 

conclude that the negative values of Aa^^ are physi- 
cally significant, and cannot be explained by conventional 
effects such as Coulomb interactions, resonance decays or 
tracking inefficiencies. 



VI. EXPERIMENT COMPARISONS 

CERN Super Proton Synchrotron (SPS) charge- 
independent <&p t measurements with a 158 GeV per nu- 
cleon Pb beam on fixed heavy ion targets (^snn ~ 17.3 
GeV) include values 0.6 ± 1.0 MeV/c for central colli- 
sions on Pb nuclei with N ~ 270 in CM pion rapid- 
ity interval 1.1 < am < 2.6 (experiment NA49) [l9| 
and 3.3 ± 0.7^'g MeV/c for central collisions on Au nu- 
clei with N ~ 162 in laboratory pseudorapidity inter- 
val 2.2 < rjiab < 2.7 (mid-pseudorapidity region) from 

the CERES experiment ^ STAR measures Act^ = 
14 ± 2 MeV/c ~ $p t for N ~ 180 when restricted to 
the CERES r\ acceptance scale [23]. All three measure- 
ments were corrected for small-scale correlations and two- 
track inefficiencies. In a following analysis |24| of the 
158 GeV per nucleon Pb-Pb fixed target collision data 
experiment NA49 reported charge independent $ Pt mea- 
surements for all charged particles in rapidity interval 
1-1 < Vn,cm < 2.6 (pion mass assumed) as a function of 
centrality. $ Pt values were found to monotonically de- 
crease from 7.2 ± 0.7 ±1.6 MeV/c for most-peripheral to 
1.4 ± 0.8 ± 1.6 MeV/c for most-central collisions. Correc- 
tions for finite two-track resolution were included, how- 
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ever the contributions of quantum and Coulomb small- 
scale correlations, estimated to be 5 ±1.5 MeV/c 0|> re ~ 

main. Quantity S pt = \J Aap[?n^ /Nftf was also reported 
by the CERES experiment |2^| with magnitude approxi- 
mately half that at STAR. Results from STAR for Ao^in 
at RHIC energy represent a striking increase over SPS 
results and markedly different centrality dependence. In 
contrast, STAR's measurement of Aa p ^^ is not signifi- 
cantly different from the NA49 result -8.5 ± 1.5 MeV/c 
in 1.1 < y n<cm < 2.6 

The PHENIX experiment at RHIC reports charge- 
independent $ Pf ss 6 ± 6 (syst) MeV/c for the upper- 
most 5% central Au-Au collision events at y/sJnv — 
130 GeV within their acceptance: < 0.35 and 

A<j) = 58.5° H3. This STAR analysis restricted to 

the PHENIX (77, 4>) acceptance scale obtained Aa Pt: n ~ 
9 ± 1 MeV/c. That value is greater than would be ex- 
pected from naive scaling from the STAR full acceptance 
scale (A?7 = 2, A<j) = 2n) to the PHENIX acceptance 
scale (Arj = 0.7, A(f> = 58.5°) [27|. The enhanced value 
for Aap^ln relative to linear scale dependence is observed 
to result from substantial nonlinear azimuth-scale (<5(/>) 
dependence of (pt) fluctuations (mainly a cos[2(</> — ^r)] 
term). 

PHENIX also reports non-zero nonstatistical (p t ) fluc- 
tuations for Au-Au collisions at ^Jsnn = 200 GeV us- 
ing quantity F Pt (2^ (proportional to <& Pt and Acr^n) 
and acceptance scales Arj = 0.7 at mid-rapidity and 
Acf> = 180° in two approximately opposed 90° spec- 
trometer arms. (p t ) fluctuations for central collisions at 
200 GeV (with two opposed spectrometer arms) are ob- 
served to be similar to those at 130 GeV (with one spec- 
trometer arm) assuming linear dependence on azimuth 
scale |23 . 

(CI CD) 

Analysis of the dependence of Aa Pt -.A on the up- 
per p t acceptance cut indicates significant contribution 
from particles with pt > 0.6 GeV/c. Subsequent studies 
of like-sign and unlike-sign two-particle correlations on 
transverse momentum space 29] for these data confirm 
that much of the observed fluctuations result from corre- 
lation excess for p t > 0.6 GeV/c. The larger magnitude 
of unlike-sign correlations relative to like-sign at higher 
p t > 0.6 GeV/c also results in Aa p ^^} < 0. These re- 
sults implicate semi-hard scattering in the initial stage 
of Au-Au collisions as a possible mechanism contributing 



to A<Tp\}^ and Aa Pt:n 
upper p t acceptance was also reported by the PHENIX 
experiment ^28| . It is therefore of interest to examine the 
predictions of available theoretical collision models which 
include hard parton scattering and/or hadronic rescatter- 
ing. 



JCD) 



Strong dependence of F Pt on the 



VII. MODEL PREDICTIONS 



HUING [5J, which incorporates p-p soft scattering and 



longitudinal color-string fragmentation phenomenology 
plus hard parton scattering and fragmentation coupled 
to a Glauber model of A-A collision geometry predicts a 

(CD 

range of Aa Pt -n up to only one- half the observed values 
in Fig. |21 hijing predictions include: 1) jet production 
enabled but without jet quenching (produces maximum 
fluctuations but still only one-half the measured values); 
2) jet production and jet quenching both enabled (vari- 
ance excess reduced by about half); and 3) no jet pro- 
duction (even smaller magnitude). In addition to un- 

derpredicting Aa p ^}n magnitudes, hijing also does not 
reproduce the observed strong centrality dependence of 
the data or the non-monotonic behavior for the more cen- 
tral collisions, but is instead approximately independent 
of centrality. 

Other collision models differ in their treatment of 
lower pt (soft) particle production, rescattering and res- 
onances, but do not include semi-hard parton scattering. 
RQMD 30 without hadronic rescattering predicts that 

(CD 

Aa Pt -n increases monotonically with centrality, reach- 
ing only half the observed value for central RHIC col- 
lisions. Initial studies of scale-dependence indicate that 
the main contribution in the RQMD model is from reso- 
nance decays and not minijets as for hijing. $ Pt pre- 
dictions from UrQMD for Au-Au collisions at y/s nn = 
200 GeV [3lJ indicate results similar to RQMD and also 
reveal strong reduction of $ Pt when hadronic rescatter- 
ing is included, rqmd and UrQMD predictions for Act^ 1 ,] 
without hadronic rescattering constitute the upper limit 
for those models. The quark-gluon string model (qgsm) 
for Pb-Pb central collisions at ^/snn — 200 GeV, when 
linearly extrapolated to the STAR acceptance scale, pre- 



dicts Act 



(CI) 
Pt-n 



10 MeV/c |32j, which is significantly less 



than the STAR measurement. 



VIII. DISCUSSION 

These fluctuation measurements, restricted to hadrons 
at lower p t (< 2 GeV/c), indicate that even central Au- 
Au collisions at RHIC are not fully equilibrated because 

(CI CD) 

A(T Pt -n would vanish for ensembles of fully equili- 
brated events (except for the relatively small contribu- 
tions from quantum and Coulomb correlations and reso- 
nance decays). Instead, Au-Au collision events at RHIC 
remain highly structured, with respect to nonstatisti- 
cal (p t ) fluctuations, as evidenced by the strong depen- 
dence on the upper p t acceptance. This result conflicts 
with assumptions underlying hydrodynamic and statis- 
tical (thermal) models conventionally applied to RHIC 
collisions. We observe no evidence of critical fluctua- 
tions associated with a possible phase transition. The 
quantity Aa p ^„ used in this analysis quantifies the sub- 
stantial differences between Au-Au collisions and simple 
models based on independent superposition of p-p colli- 
sions. We have demonstrated that the observed charge- 
independent and charge-dependent nonstatistical fluctu- 
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ations cannot be explained in terms of final-state quan- 
tum and Coulomb correlations and resonance decays or in 
terms of experimental effects such as two-track inefficien- 
cies and time dependences of experimental apparatus. 

( CI) 

The observed strong energy dependence of Acr pt: „ from 
SPS to RHIC and the failure of conventional theoretical 
models to describe these new RHIC fluctuation data in- 
dicate that significant new dynamical mechanisms play 
a role in Au-Au collisions at RHIC, mechanisms which 
substantially affect the correlation structure of final-state 

fen 

transverse momentum. The increase of Aerp t: „ with p t 
upper limit, combined with apparent saturation and even 
reduction of Aa^ r l for the more central Au-Au colli- 
sions, suggests that semi-hard parton scattering and sub- 
sequent dissipation of parton momentum by coupling to 
an increasingly dense, possibly colored medium may ac- 
count for these observations. Detailed studies of correla- 
tion structure in both transverse and longitudinal mo- 
mentum components will be reported in the near fu- 
ture 0,111113. 

IX. CONCLUSIONS 

This first large-acceptance measurement of (pt) fluctu- 
ations at RHIC reveals intriguing deviations from a sta- 
tistical reference. We observe a 13.7 ± 1.4% (stat+syst) 
r.m.s. fractional excess of charge-independent fluctua- 
tions in \^n((pt) — pt) (17 ± 2% (stat+syst) when ex- 
trapolated to 100% of primary charged hadrons in the 
STAR acceptance) for the 15% most-central events which 
varies smoothly and non-monotonically with centrality. 
This observation of strong nonstatistical (p t ) fluctua- 
tions demonstrates that RHIC events are not fully equi- 
librated, even in the lower p t sector for central events, 
contradicting a basic assumption of hydrodynamic and 
statistical models. There is no significant evidence for 
anomalous event classes as might be expected from criti- 
cal fluctuations. Comparisons with SPS experiments in- 
dicate that charge-independent fluctuations are qualita- 
tively larger at RHIC, whereas charge-dependent fluctu- 
ations are not. A PHENIX result at 130 GeV for charge- 
independent fluctuations, compatible with zero within 
their systematic error, is consistent with a significant 
non-zero STAR measurement restricted to the PHENIX 
acceptance. Based upon studies of the higher p t contri- 
bution and various model comparisons we speculate that 
these fluctuations may be a consequence of semi-hard 
initial-state scattering (minijets) followed by parton cas- 
cading in the early stage of the Au-Au collision which 
is not fully equilibrated prior to kinetic decoupling |34j . 
Such strong fluctuations have not been observed previ- 
ously in heavy ion collisions and are at present unex- 
plained by theory, thus pointing to the possibility of new, 
or perhaps unexpected dynamical processes occuring at 
RHIC. Identification of the dynamical source(s) of these 
nonstatistical fluctuations is underway p9l | and will con- 
tinue to be studied in the future. 
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APPENDIX A: 

In this appendix the total variance is defined. The 
scale invariance of total variance, an alternative state- 
ment of the central limit theorem [l5|, [l6j , then motivates 
the definition of fluctuation measure Aap[^ used in this 
analysis. 

A detector acceptance (Ar7,A</>) (generically Ax) on 
axial momentum space (77, 4>) can be divided into bins 
of size (5r],5(j)) (generically Sx). Each bin then contains 
event- wise scalar p t sum 

n a (Sx) 

p t , a (Sx) = ^ Pt,ah ( Al ) 
i-1 

where a is a bin index and n a (Sx) is the event-wise mul- 
tiplicity in bin a. Fluctuations in p t ^ a (Sx) relative to 
n a {Sx) could be measured by the variance of the ratio 
(pt)a = Pt,a(Sx) /n a {5x). However, to minimize contri- 
butions from event-wise and bin-wise variations in n a (Sx) 
(a source of systematic error) we instead compute the to- 
tal variance of difference pt, a {Sx) — n a (Sx)pt, defined by 

M(Ax,Sx) 

Y, 2 pt:n (Ax,6x)= (PtA^)-n a (Sx)pt) 2 ,(A2) 

a=l 

where M(Ax,Sx) is the event-wise number of occupied 
bins of size Sx in acceptance Ax and the bar denotes an 
average over all events. Typically M(Ax, 5x) = Ax/Sx, 
a constant for all events except when 5x <C Ax and some 
bins are unoccupied. 

For the analysis described in this paper we are inter- 
ested in two limits of Eq. (|A2() . the acceptance scale 
5x = Ax with M — 1, and a 'single-particle' scale 
Sx <C Ax such that each occupied (77, 4>) bin contains 
a single particle, with M — > n(Ax) = N(Ax), the event- 
wise total multiplicity in the acceptance. For a collection 
of reference events (c/ Sec. IIII(I obtained by independent 
Pt sampling from a fixed parent distribution (also referred 
to here as 'CLT conditions') quantity T,p t . n (Ax, Sx) is in- 
dependent of bin size Sx. We illustrate this scale invari- 
ance under CLT conditions for the above two limits and 
for arbitrary scale Sx as follows. 
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In the single-particle scale limit each occupied bin con- 
tains only one particle, and the bin index is equivalent to 
a particle index: pt. a (Sx) — > pt t i (transverse momentum 
of particle i) and n a {8x) — » 1. £p t . n (Ax, Sx) then has 
the limit 



Y; 2 pt . n {Ax,6x<t: Ax) iV(Ax)a| t , 



(A3) 



where N(Ax) is the mean total event multiplicity, and 
the variance of the inclusive pt distribution is explicitly 

4 = EiJt X) (Pt,i - Pt) 2 /N(Ax). In the limit Sx -» Ax, 
M(Ax, 5x) — > 1, the event-wise single-bin occupancy is 
N(Ax), and T, 2 . n (Ax, Sx) becomes 



T 2 



: „(Ax, Sx = Ax) 



PnN{A. 

N(Ax) 



where the sum includes all values of event multiplicity 
7V(Ax) represented in the event ensemble, pn = Eat/e is 
the fraction of events in the ensemble with multiplicity 
iV(Ax), and ct? > = ((pt)N — Pt) 2 is the variance of the 
(pt) distribution for the subset of events with multiplicity 
iV(Ax). If CLT conditions apply, then 

X 2 pttn (Ax,Ax) (C = T) J2PNN(Ax)a 2 t (A5) 

JV 

= N(Ax)a 2 pt , 

where CLT relation aj > = c? IN (Ax) was invoked. 

The equivalence under CLT conditions of Sp 4 . n (Ax, Sx) 
for these two limiting scale values is thus established. 

Generalizing the latter argument, the total variance at 
arbitrary scale Sx in Eq. I|A2() can be reexpressed as 



M(Ax,8x) 

E 2 (Ax,Sx) = V nl(5x)((p t ) a -p t ) 2 (A6) 



a=l 



= M(Ax,Sx) p n n 2 (Sx)<j 2 {pt)n , 

n{8x) 

where sums over events and bins were rearranged as sums 
over bin-wise multiplicity n(Sx) and over bins a which 
have that value of multiplicity, p n is the fraction of bins 
in the event ensemble with multiplicity n(Sx), and P 2 pt ) 
is the variance of (p t ) — pt within that subset of bins 



r <P 



t ) n = ((Pt)u-Ptf 



(A7) 



The bar in Eq. (|A7|) indicates an average over all bins 
in the event ensemble with multplicity n. For CLT 



Op t /n(5x) for any n, and, since 



conditions cr? , 

(Pt)r, 

M(Ax,5x)n(Sx) = N(Ax), Eq. (|A6l) therefore becomes 



E 2 pt .JA Xl Sx) = N(Ax)al, 



(A8) 



which demonstrates the general scale invariance of 
Y>l t , n (Ax,5x) for CLT conditions. 

Deviations from central limit conditions signal the 
presence of two-particle correlations (e.g., p t samples are 
not independent). The total variance is then no longer 
scale invariant, and its scale dependence reflects the de- 
tailed structure of those correlations. We therefore define 
a total variance difference between arbitrary scales 5x\ 
and 8x2, where 8x\ < 8x2, as 



AY, 2 pt . n (Ax,5x 1 ,5x 2 ) 



= Kf.n( Ax ^ X 2)- 



t (Ax,Sxi 



(A9) 



where AE?_. n (Ax, Sxi, 8x2) = if CLT conditions apply 
in the scale interval [Sx\, 8x2]- 

The total variance difference depends by construction 
on the detector acceptance (and on the collision system or 
particpant number). We can remove those dependences 
in several ways whose choice depends on the physical 
mechanisms producing the correlations. For this applica- 
tion we divide by the total multiplicity in the acceptance 
to obtain a fluctuation measure per final-state particle. 

If CLT conditions are approximately valid, n(Sx)a 2 pt ^ 
in Eq. (|A6|) is nearly constant and can be removed from 
the weighted summation over n, resulting in 



E^JAx.Sx) ~ N(Ax)n(Sx)(( Pt )-p t ) 2 , (A10) 

a factorized form in which acceptance and scale depen- 
dences are separated. The total variance difference for 
8x2 — Sx and Sx\ <C Ax is then given by 

AEp t: „(Ax,<5xi < Ax,8x) 



N(Ax) x [n(8x)((p t ) - p t ) 2 
N(Ax)Ao 2 (5x) 



(All) 



combining Eqs. (|A3|I and ljA10|l . In Eqs. l|A10f> and (|A11|) 

the bar denotes an event-wise sum over occupied bins and 
an average over all events. The (pt) fluctuation excess 

measure Ao-pl^ in Eq.® is thus identified as the total 
variance difference in Eq. (|A11I) per final-state particle, 
evaluated at the acceptance scale Sx = Ax. 
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